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We show how to simplify the computation of the entangle- 
ment of formation and the relative entropy of entanglement 
for states, which are invariant under a group of local symme- 
tries. For several examples of groups we characterize the state 
spaces, which are invariant under these groups. For specific 
examples we calculate the entanglement measures. In par- 
ticular, we derive an explicit formula for the entanglement of 
formation for U ® {/-invariant states, and we find a counterex- 
ample to the additivity conjecture for the relative entropy of 
entanglement. 



I. INTRODUCTION 

One of the reasons the general theory of entangle- 
ment has proved to be so difficult is the rapid growth 
of dimension of the state spaces. For bipartite entan- 
glement between d\- and c?2-dimensional Hilbert spaces, 
entanglement is a geometric structure in the [<i\d\ — 1)- 
dimensional state space. Hence even in the simplest non- 
trivial case (di = c?2 = 2; 15 dimensions) naive geometric 
intuitions can be misleading. On the other hand, the 
rapid growth of dimensions is partly responsible for the 
potential of quantum computing. Hence exploring this 
complexity is an important challenge for quantum infor- 
mation theory. 

Model studies have been an important tool for devel- 
oping and testing new concepts and relations in entangle- 
ment theory, both qualitative and quantitative. In this 
paper we explore a method for arriving at a large class 
of models, which are at the same simple, and yet show 
some of the interesting features of the full structure. 

The basic idea, namely looking at sets of states which 
are invariant under a group of local unitaries is not new, 
and goes back to the first studies of entanglement in 
the modern sense. Two classes, in particular, have been 
considered frequently: the so-called Werner states (after 
M), which are invariant under all unitaries of the form 
U (g) U, and the so-called isotropic states J3|, which are 
invariant under all U®U, where U is the complex conju- 
gate of U in some basis. Symmetry has also been used in 



*Electronic Mail: k.vollbrecht@tu-bs.de 
^Electronic Mail: R.Werner@tu-bs.de 



this way to study tripartite entanglement ||. A re- 
cent paper of Rains || discusses distillible entanglement 
under symmetry, so we have eliminated the pertinent re- 
marks from this paper. 

Several of the ingredients of our general theory, for ex- 
ample the role of the twirl projection and the commutant, 
have been noted in these special cases and can be consid- 
ered to be well-known. The computation of the relative 
entropy of entanglement was known Js) for Werner 
states. The first study in which symmetry is exploited to 
compute the entanglement of formation || beyond the 
Wootters formula is |Q , where the case of isotropic 
states is investigated. Our theory of entanglement of for- 
mation can be viewed as an abstract version of arguments 
from that paper. 

What is new in the present paper is firstly the gener- 
ality. We regard our theory as a toolkit for constructing 
examples adapted to specific problems, and we have tried 
to present it in a self-contained way facilitating such ap- 
plications. Exploring all the possibilities would have been 
too much for a single paper but, of course, we also have 
some new results in specific examples. 

Our most striking specific result is perhaps a coun- 
terexample to the conjecture that the relative entropy 
of formation should be additive. The evidence in favor 
of this conjecture had been partly numerical, but it was 
perhaps clear that a random search for counterexamples 
was not very strong evidence to begin with: the relative 
entropy of entanglement is defined by a variational for- 
mula in a very high dimensional space, whose solution is 
itself not easy to do reliably. In addition, the additivity 
conjecture is true on a large set in the state space, so 
unless one has a specific idea where to look, a random 
search may well produce misleading evidence. The sec- 
ond strong point in favor of the additivity conjecture had 
been a Theorem by Rains (Theorems 4 and 5 in |H| ) im- 
plying a host of non-trivial additivity statements. How- 
ever, our counterexample satisfies the assumptions of the 
Rains' Theorem, so that Theorem is, unfortunately, false. 

Further specific results in our paper are the formulas 
for entanglement of formation and relative entropy of en- 
tanglement for Werner states. 

The paper is organized as follows: In Section II we re- 
view the essential techniques for the investigation of sym- 
metric states and describe how the partial transposition 
fit in this context. Section II. D presents a zoo of different 
symmetry groups. Some of these are used later, others 
are only presented as briefly, to illustrate special prop- 
erties possible in this setup. We hope that this list will 
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prove useful for finding the right tradeoff between high 
symmetry, making an example manageable, and richness 
of the symmetric state space, which may be needed to see 
the phenonmenon under investigation. In Section ^ we 
briefly recapitulate the definitions of the entanglement of 
formation and the relative entropy of entanglement and 
the additivity problem. In Section IV we turn to the 
entanglement of formation. We show first how the com- 
putation may be simplified using local symmetry. These 
ideas are then applied to the basic symmetry groups UU 
and UU, arriving at an explicit formula in both cases 
(the results for UU are merely cited here for complete- 
ness from work of the first author with B. Terhal [p2|). 
For the group OO of orthogonal symmetries, which uni- 
fies and extends these two examples, we find formulas in 
large sections of the state space. Section V deals with 
the relative entropy of entanglement. Again we begin by 
showing how the computation is simplified under sym- 
metry. We then present the counterexample to additivity 
mentioned in the introduction. Some possible extensions 
are mentioned in the concluding remarks. 



II. SYMMETRIES AND PARTIAL TRANSPOSES 

From the beginning of the theory of entanglement the 
study of special subclasses of symmetric states has played 
an important role. In this section we give a unified treat- 
ment of the mathematical structure underlying all these 
studies. For simplicity we restrict attention to the bipar- 
tite finite dimensional case, although some of the gener- 
alizations to more than two subsystems Q and infinite 
dimension are straightforward. So throughout we will 
consider a composite quantum system with Hilbert space 
H = Hi <g> H.2, with dim Hi = d% < oo. We denote the 
space of states (=density operators) on H as S(H), or 
simply by S. The space of all separable states (explained 
in subsection II B) is denoted as T>. 



A. Local symmetry groups 

Two states p, p' are regarded as "equally entangled" 
if they differ only by a choice of basis in Hi and Hi or, 
equivalently, if there are unitary operators U acting onHi 
such that p' = (Ui ® U-i)p{U\ ® U2)* ■ If in this equation 
p' = p, we call U — (Ui ® U 2 ) a (local) symmetry of the 
entangled state p. Clearly, the set of symmetries forms a 
closed group of unitary operators on Hi ® Hi- We will 
now turn this around, i.e., we fix the symmetry group 
and study the set of states left invariant by it. 

So from now on, let G be a closed group of unitary 
operators of the form U = (U\ ® U%). As a closed sub- 
group of the unitary group, G is compact, hence car- 
ries a unique measure which is normalized and invariant 
under right and left group translation. Integrals with 
respect to this Haar measure will just be denoted by 



"/ dU " , and should be considered as averages over the 
group. In particular, when G is a finite group, we have 
JdU f(U) = |Gr 1 Ec/ eG /( [/ )- An important ingredi- 
ent of our theory is the projection 
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for any operator A on Hi ® H2, which is sometimes re- 
ferred to as the twirl operation. It is a completely pos- 
itive operator, and is doubly stochastic in the sense that 
it takes density operators to density operators and the 
identity operator to itself. Using the invariance of the 
Haar measure it is immediately clear that "PA = A" is 
equivalent to U [U, A] = for all U G G" . The set of all 
A with this property is called the commutant of G. We 
will denote it by G' , which is the standard notation for 
commutants in the theory of von Neumann algebras. It 
will be important later on that G' is always an algebra 
(closed under the operator product), although in general 
P(AB) ^ (PA)(PB). Computing the commutant is al- 
ways the first step in applying our theory. Typically, one 
tries to pick a large symmetry group G from the out- 
set, so the commutant becomes a low dimensional space, 
spanned by just a few operators. 

Our main interest does not lie in the set G 1 of G- 
invariant observables, but dually, in the G-invariant den- 
sity operators p with Pp — p. As for observables this set 
is the projection PS of the full state space under twirling. 
The relation between invariant observables and states is 
contained in the equation 



tx{P{p)A) = tr(pP(A)) , 



(2) 



which follows easily by substituting U 1— ► U* = U^ 1 
in the integral ([!]), and moving one factor U under the 
trace. Due to this equation, we do not need to know 
the expectations tr(pA) for all observables A in order 
to characterize a G-invariant p, but only for the invari- 
ant elements P(A) £ G' . Indeed, if we have a linear 
functional / : G' — ► C, which is positive on positive op- 
erators, and normalized to /(I) = I , that is a state on 
the algebra G 1 in C*-algebraic terminology, the equation 
tr(pA) = f(P(A)) uniquely defines a G-invariant density 
operator p, because P preserves positivity and P(l) = 1. 
Under this identification of G-invariant density operators 
and states on G' it becomes easy to compute the image 
of a general density operator p under twirling. Using 
again equation (||) we find that Pp is determined simply 
by computing its expectation values for A <E G' , i.e., its 
restriction to G'. 

Let us demonstrate this in the two basic examples of 
twirling. 

Example 1: The group UU (Werner states). 

We take the Hilbert spaces of Alice and Bob to be the 
same (H = Hi ® Hi), and choose for G the group of all 
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unitaries of the form U ®U , where U is a unitary on Hi- 
As an abstract topological group this is the same as the 
unitary group on Hi, so the Haar measure on G is just 
invariant integration with respect to U. It is a well-known 
result of group representation theory, going back to Weyl 
fl4| or further, that the commutant of G is spanned by 
the permutation operators of the factors, in this case the 
identity 1 and the flip defined by F((p eg) ip) = ip eg) (p, or 
in a basis \i) of Hi, with \ij) = \i) ® \j): 



F 



(3) 



h3 



Hence the algebra G' consists of all operators of the form 
A = al + f3F. As an abstract *-algebra with iden- 
tity it is characterized by the relations F 2 — 1 and 
F* = F. Thus G-invariant states are given in terms 
of the single parameter tr(pF), which ranges from —1 to 
1. Note that an invariant density operator can be written 
as p = al + bF with suitable a,i £ R. But as we will 
see, the parameters a, b are less natural to use, and more 
dimension dependent than tr(pF). 

Example 2: The group UU (isotropic states). 

Again we take both Hilbert spaces to be the same, and 
moreover, we fix some basis in this space. The group G 
now consists of all unitaries of the form U®U, where U is 
a unitary on Hi, and U denotes the matrix element- wise 
complex conjugate of U with respect to the chosen basis. 
One readily checks that the maximally entangled vector 
$ = Y^i I**) is invariant under such unitaries, and indeed 
the commutant is now spanned by 1 and the rank one 
operator 



F = |3X*l=X>>0'il- 



(4) 



This operator is positive with norm d = ||$|| = dim Hi, 
so the invariant states are parametrized by the interval 
[0, d\. These claims can be obtained from the first exam- 
ple by the meth od of partial transposition discussed in 
Subsection II C. 



It is perhaps helpful to note that there are not so many 
functions U > U , taking unitaries on Hi to unitaries on 
the same space Hi, such that the operators of the form 
U <g)U again form a group. For this it is necessary that 
U i i U is a homomorphism, so, for example, U = U* 
does not work. Inner homomorphisms, i.e., those of the 
form U = VUV* are equivalent to Example 1 by a trivial 
basis change in the second factor, given by V. Similarly, 
functions differing only by a scalar phase factor give the 
same transformations on operators, and should thus be 
considered equivalent. Then (up to base changes and 
phase factors) all functions U i— > U not equivalent to the 
identity are equivalent to Example 2, i.e., the above list 
is in some sense complete. However, many interesting 



examples arise, when the Hilbert spaces are not of the 
same dimension, or the group of operators in the first 
factor is not the full unitary group. 

Computing PS in Examples 1 and 2 is very simple, 
because it is just an interval. We will encounter more 
complicated cases below, in most of which, however, the 
algebra G' is abelian. When G' has dimension k, say, it 
is then generated by k minimal projections, which corre- 
spond precisely to the extreme points of PS. Therefore 
the state space is a simplex (generalized tetrahedron). 



B. How to compute the separable states PX> 

For the study of entanglement of symmetric states it 
is fundamental to know which of the states in PS are 
separable or "classically correlated" Q|, i.e., convex com- 
binations 



P 



(5) 



of product density operators. We denote this set of states 
by V. Because we assume the group G to consist of 
local unitaries, it is clear that for a separable state p 
the integrand of Pp consists entirely of separable states, 
hence Pp is separable. Hence PT> C (T> (1 PS). But here 
we even have equality, because any state in (T> (~l PS) is 
its own projection. Hence 



V n PS = PV 



(6) 



In order to determine this set, recall that by decomposing 
Pi*2 m © P ure states, we may even assume the 

Pi°2 in (^) to be pure. If we compute Pp termwise, we 
find that each p € PT> is a convex combination of states 
P(<7i (8><72) with pure cr, = \(j>i)(4>i\. Thus we can compute 
PV in two stages: 

• Choose a basis in G' , consisting, say of k hermi- 
tian operators A a and compute the expectations of 
these operators in arbitrary pure product states: 

a n = (<t>-\ ® <t>i\A n \6-\ ® 6o) , 



this determines the projections P(ai (£> 02). 

• Determine the set of real fc-tuples (01, . . . , a^) ob- 
tained in this way, as the range over all normal- 
ized vectors. 

• Compute the convex hull of this set. 

Two simplifications can be made in this procedure: 
firstly, we always have I e G' , so by choosing Ak = 1, 
it suffices to work with the (k — l)-tuples (01, . . . , cifc-i). 
Secondly, the vectors (pi <g> (pi and U((pi ® fa) with U S G 
give the same expectations, so when determining the 
range one can make special choices, as long as one vector 
is chosen from each orbit of product vectors under G. 
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Let us illustrate this procedure in the two basic exam- 
ples above: In Example 1 we only need to compute 

(0®</>|F|<^} = l«#)| 2 - (7) 

Clearly, this quantity ranges over the interval [0, 1], and 
a [/[/-invariant state p is separable iff tr(pF) > |Q. 
Similarly, in Example 2: 

i 

which again ranges over the interval [0, 1] . Note, however, 
that the state space in this case is the interval [0, d}. The 
fact that the two state space intervals [—1,1] for UU and 
[0, d] for UU intersect precisely in the separable subset 
[0, 1] is an instance of the Peres-Horodecki criterion for 
separablility, as we now proceed to show. 

C. Partial transposition 

The partial transpose of an operator on Hi ® Ji 2 is 
defined in a product basis by transposing only the indices 
belonging to the basis of 7i 2 , and not those pertaining 
Tlx- Equivalently, we can define this operation as 

e 2 {A(g>B) = A(g> 0(B) , (9) 

where 6(1?) denotes the ordinary matrix transpose of 
B. This also depends on the choice of basis in Tl 2 , so 
from now on we assume a basis of Tl 2 to be fixed. This 
equation suffices to define O2, because all operators on 
H-i <8> Ti-2 can be expanded in terms of product operators. 
The partial transpose operation has become a standard 
tool in entanglement theory with the realization that the 
partial transpose of a separable density operator is again 
positive. This is evident from Equations (|5|) and (|^), and 
the observation that the transpose of a positive operator 
is positive. In 2 £g> 2 and 2 ® 3 Hilbert space dimensions, 
this criterion, known as the Peres-Horodecki criterion, 
is even sufficient for separability [ jl5| . For all higher di- 
mensions sufficiency fails in general. States with posi- 
tive partial transpose ( "ppt-states" ) are known not to be 
distillible, i.e., even when many copies of such a state 
are provided, it is not possible to extract any highly en- 
tangled states by local quantum operations and classical 
communication alone. 

For special classes of states on higher dimensional 
Hilbert spaces the ppt-property may still be sufficient for 
separability. Pure states are a case in point, and so are 
some of the spaces of symmetric states studied in this pa- 
per. Let us check how the action of a product unitary is 
modified by partial transposition. If U, Ai are operators 
on Hi (i = 1, 2), we find 

e 2 (([/i ® U 2 )( Ai (8 A 2 ){Ul ® E/J)) 

-e 2 (([/ 1 A 1 [/*)® (L/ 2 A 2 [/*)) 



= (iTiAitf?) ® (e([/ 2 *)e(A 2 )e([/ 2 )) 

- {U 1 ®U 2 )®2{Ai® A 2 ){U X ® U 2 )* . 

Note that by linearity we can replace in this equation 
A\ <g> A 2 by any other operator on Tlx % Ti 2 . This compu- 
tation motivates the following definition: For any group 
G of product unitaries we denote by G the group of uni- 
taries U\ ®U 2 , where U\®U 2 G G. For example, for 
G = UU of Example 1 we get G = UU, and conversely. 

There is a slightly tricky point in this definition, be- 
cause the map U\ ® U 2 1— > Ux ® U 2 is not well defined: If 
we multiply U\ by a phase and U 2 with the inverse phase, 
the operator U\ ® U 2 does not change, but U\ ® U 2 picks 
up twice the phase. What the definition therefore re- 
quires is to take in G all operators arising in this way. 
Repeating the "twiddle" operation may thus fail to lead 
back to G, but instead leads to G enlarged by the group 
of phases. It is therefore convenient to assume that all 
groups under consideration contain the group of phases. 
We may do so without loss of generality, since the phases 
act trivially on operators anyhow, and hence the twirling 
projection P is unchanged. 

If we integrate the above computation with respect to 
a group G of local unitaries, and introduce P for the 
twirling projection associated with G, we get the funda- 
mental relation 

e 2 p = pe 2 . (10) 

Since 6 2 is a linear bijection on the space of all operators 
on Tlx <g> Tl 2 , we immediately find the relations between 
the ranges of P and P: 

e 2 (C) = g' , (11) 

i.e., the operators invariant under G are precisely the 
partial transposes of those invariant under G. This has 
a surprising consequence: taking the partial transposes 
of an algebra of operators in general has little chance 
of producing again an algebra of operators, since 9 2 is 
definitely not a homomorphism. That is, in general one 
would not expect that the operator product of two partial 
transposes is again the partial transpose of an element of 
the original algebra. If the algebra arises as the commu- 
tant of a group of local unitaries, however, we get again 
a commutant, hence an algebra. 

The first application of Equation (|ll|) is the computa- 
tion of the commutant in Example 2: With G = UU we 
find the partial transposes of the operators in G', i.e., the 
operators 9 2 (al + (3F) = al + 0, since 62(F) = F. 

Another application is the determination of the set of 
ppt-states. Cne might think that a special form for p, 
entailed by its G-invariance, is not necessarily helpful for 
getting spectral information about 2 p. However, since 
6 2 G" is an algebra, and often enough an abelian one, 6 2 p 
is, in fact, easily diagonalized. 

A good way to represent this connection is to draw the 
state spaces of G and G (i.e., PS and PS) in the same 
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diagram. Since in general G' and G' need not intersect 
except in the multiples of the identity (see Examples 1 
and 2), the projected state spaces PS and PS in gen- 
eral have only the trace-state in common. Hence they 
don't fit naturally in the same diagram. However, the 
partial transposes of PS lie in G', more precisely in the 
hyperplane of hermitian elements with trace 1. The same 
hyperplane contains PS. In the pair of Examples 1 and 
2, we get Figure 1. 



(F) 



PS 



PV 



PV 



PS 



(F) 



-1 -0.5 0.5 1 d 

FIG. 1. The state spaces for Werner and isotropic states 
are just intervals. Drawn together in one diagram the inter- 
section gives us the space of PPT-states, which is in this case 
equivalent to the separable space. 



Note that by exchanging the roles of G and G, we 
get exactly the same diagram, up to maybe an affinc 
transformation due to a different choice of coordinates: 
the two diagrams are simply related by taking partial 
transposes. When G and G are swapped in this way, 
the picture of PT> remains correct: since Q2D — T>, it 
suffices to compute the projection of the separable subset 
for G. By definition, the intersection of PS and Q2PS 
is the convex set of G-invariant ppt-states. It always 
contains PT>, but this inclusion may be strict. In the 
simple case of Figure 1 PV = PS n Q2PS, which is the 
same as saying that the Peres- Horodecki criterion is valid 
for states invariant under either G or G. 



D. Further examples of symmetry groups 



Example 3: Orthogonal groups: G = OO. 

The two basic examples can be combined into one by 
taking the intersection of the two groups: G = UU H 
UU this is the same as the subgroup of unitaries U ®U 
such that U — U, i.e., such that U is a real orthogonal 
matrix. Clearly, both the [/[/-invariant states and the 
U [/-invariant states will be G-invariant, so we know that 
G' is at least^the algebra generated by UU' and UU , 
i.e., it contains 1,F, and F. Since FF = FF = F, the 



linear span of these three is already an algebra, and is 
spanned by the minimal projections 



Po 



-F 



Pl = -(1-F) 

P2 = -(i-F) - If 



(12) 
(13) 
(14) 



which corresponds precisely to the decomposition of a 
general 3 x 3-matrix into multiple of the identity, an- 
tisymmetric part, and symmetric traceless part. This 
decomposition of tensor operators with respect to the or- 
thogonal group is well known, so we have identified G' . 

The extremal G-invariant states corresponding to these 
three minimal projections are plotted in Figure @ in a co- 
ordinate system whose axes represent the expectations 
of F and F, respectively. The plane of this drawing 
should be considered as the hermitian G-invariant oper- 
ators of trace one. This plane is mapped into itself by 
partial transposition (since G = G) , and the coordinates 
are chosen such that partial transposition is simply the 
reflection along the main diagonal. 



(F) 




(F) 



FIG. 2. State spaces for OO and OO invariant states plot- 
ted for d = 3. The UU and UU invariant states are drawn as 
thin lines. 



1 In general, the commutant (GClH)' may be properly larger 
than the algebra G' V H' generated by G' and H' . The equa- 
tion (.An B)' = A' VB' is valid only for algebras, and follows 
readily from the equation (A' V B')' = (A" n B"), and the 
bicommutant theorem JTt| , which characterizes M" as the 
algebra generated by M. However, the algebras G" and H" 



may have an intersection, which is properly larger than the al- 
gebra generated by their intersection. For example, for any ir- 
reducible represented group G" is the algebra of all operators, 
but two such groups may intersect just in the identity. Hence 
some caution has to be exercised when computing (G n H)' 
for general groups. 
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The intersection of PS and P6> is the square [0, 1] x 
[0,1]. Is the Peres- Horodecki criterion valid for these 
states? All we have to do to check this is to try to get 
some pure product states, whose expectations of F and 
F fall on the corners of this square. For a product vector 
cj) ® ip we get the pair of expectations 

(|(#/>)l 2 , MW) ■ 

Here ip denotes the complex conjugate of ip in a basis in 
which the representation is real. Now the point (1, 1) in 
the square is obtained, whenever <p = ip is real, the point 
(0, 0) is obtained when <p and ip are real and orthogonal, 
and the point (1, 0) is obtained when ip — <p, and (<p\<p) — 
0, for example <fi — (1, *, 0)/ \/2. Symmetrically we get 
(0,1) with the same <j> and ip — <j>. Hence all four corners 
are in PT>, and as this is a convex set we must have 

PV = PSne 2 (PS). 

Example 4: SU2-representations. 

A class of examples, in which arbitrary dimensions of Hi 
and Ti.2 can occur is the following. Let u i— > T> J U denote 
the spin j irreducible representation of SU 2 - Then we 
can take 

G= {Di 1 <E>V^\ue SU 2 } , (15) 

where (2jk + 1) is the dimension of Hk (k = 1,2). Since 
the jk also take half-integer values, these dimensions can 
be any natural number > 1. It is known from just 
about any quantum mechanics course (under the key 
word "addition of angular momenta") that the tensor 
product representation 2? J1 g) I?- 72 is decomposed into 
the direct sum of the irreducible representations V s with 
s = \ji ~ h\, \ji -h\ + 1 ! • • • (ii +h), each of these repre- 
sentations appearing with multiplicity 1. Therefore, the 
commutant of G is spanned by the projections onto these 
subspaces, and is an abelian algebra. 

Note that since the spin-1 representation of SU2 is the 
orthogonal group in 3 dimensions, the case j\ = j 2 = 1 
corresponds precisely to the previous example with d = 3. 
We have no general expression for the separable subsets, 
nor even for the partially transposed sets in these exam- 
ples. We believe, however, that this class of examples 
deserves further investigation. 

Example 5: Bell diagonal states. 

In this example we show that the group G can also be 
abelian, and we make contact with a well investigated 
structure of the two qubit system. So let Hi = H 2 = C 2 , 
and let <7k, k = 1, 2, 3 be the Pauli matrices, and <jq = I. 
Then the set 

G = {I, -<7i <8> 0-1, -02® 02, -0-3 (g) cr 3 } (16) 

forms a group, which is isomorphic to the Klein 4-group, 
and abelian (G C G 1 ). It is even maximally abelian, 



i.e, the algebra G" generated by G is equal to, and not 
just contained in G' . The minimal projections in G' are 
|*fcX**li fc = 0, 1, 2, 3, where the * fc are the magical Bell 
Basis |@: * = (|ll) + |22))/v/2, and * fc = z(l®a fc )* 
for fc = 1, 2, 3. In this basis the group elements and their 
negatives are the diagonal operators with diagonal ele- 
ments ±1, of which an even number arc — 1. Hence the 
G-invariant states are the tetrahedron of density opera- 
tors which are diagonal in Bell Basis. 

The partial transpose is easy to compute: only <j 2 
changes sign under transposition. Hence if we draw the 
state space in a coordinate system, whose three axes 
are the expectations of the group elements — Ok ® &k 
(k = 1,2,3), the Bell states are the corners (1,1,1), 
(1, —1, —1), (—1, —1, 1), and (—1, 1, —1) of the unit cube, 
from which their partial transposes are obtained by mir- 
ror reflection X2 1— > — X2- That is, the partially trans- 
posed states occupy the remaining four corners of the 
unit cube. The ppt-subset, which is equal to the separa- 
ble subset since we are in 2 g) 2 dimensions, is hence the 
intersection of two tetrahedra, and is easily seen to be an 
octahedron. 




FIG. 3. State spaces for Bell diagonal states. 



Example 6: Finite Weyl Systems. 

In the examples so far the groups G and G were iso- 
morphic or even equal. In this example, which extends 
the previous one, we see that the two groups and their 
commutants can be very different. 

We let d be an integer, and introduce on C d the Weyl 
operators, given by 

W(x,y)\z)=u; xz \z-y) , (17) 

where to = exp(27ri/p). These are unitary, and satisfy 
the 'Weyl relations" 



G 



W(xuyi)W(x2,m) =u- x ^W( Xl +x2,yi + y2) • (18) 

Hence these operators, together with the p th roots of 
unity form a group. On <C d ® C d we introduce the op- 
erators W(x 1 ,y 1 ,x 2 ,y 2 ) = W(x 1 ,y 1 ) ®W{x 2l y 2 ), and 
take 

G = {uj z W(x,y,x,y) \x,y,z = 0, . . . , d - 1} . (19) 

The commutant is readily computed from the Weyl rela- 
tions to be 



G' = span{vy(x, y, -x, -y) \x, y = 0, . 



1} ■ (20) 



The Weyl operators in G' satisfy Weyl relations with u> 
replaced by to 2 . If cZ is odd, such relations are equivalent 
to the Weyl relations (^8[) for a d-dimensional system, 
and hence G' is isomorphic to the d x d-matrices. 

On the other hand, complex conjugation of W(x,y) 
just inverts the sign of x, so G contains the Weyl oper- 
ators W(x,y,—x,y). But this time, rather than getting 
twice the Weyl phase, the phases cancel, and G is abelian. 
One also verifies that 

G' = span{ W(x, y, -x, y) \x, y = 0, . . . , d - l} . (21) 

is spanned by G, so this algebra is even maximally 
abelian: it contains d 2 one-dimensional projections, 
which thus form the extreme points of PS. Hence we get 
the following picture: the set PS of G-invariant states is 
isomorphic to the space of d x d-density operators, and 
the G-invariant operators with positive partial transpose 
are a simplex spanned by 9 extreme points, which are 
mapped into each other by the action of a d x d Weyl 
system. The intersection PS n 82 (PS) is a rather com- 
plicated object. We do not know yet whether it differs 
from PT>. 

Example 7: Tensor products. 



Additivity problems for entanglement (sec Section III C 
for a brief survey) concern tensor products of bipartite 
states, which are taken in such a way as to preserve the 
splitting between Alice and Bob. Thus in the simplest 
case we have four subsystems, described in Hilbert spaces 
Hi,/Ci, i — 1,2, such that systems Tli and ICi belong 
to Alice, systems H 2 and JC 2 belong to Bob, and such 
that the systems in Hi are prepared together according 
to a density matrix p on Tii C§> Tl 2 and, similarly, the 
remaining systems are prepared according to a, a density 
operator on JCi^JC 2 . We wish to study the entanglement 
properties of p®a, when both these density matrices are 
assumed to be invariant under suitable groups of local 
unitaries. 

Let us denote by G (resp. H) the group of local uni- 
taries on "H\®'H 2 (resp. by K.\®K. 2 ), and assume p and a 
to be invariant under the respective group. Then, clearly, 
p ® a is invariant under all unitaries U\ ® U 2 <8> V\ ® V 2 , 



where U\®U 2 £ G and V% ® V 2 £ H. These again form a 
group of local unitaries, denoted by G®H, where "local" 
is understood in the sense of the Alice— Bob splitting of 
the system, i.e., the unitary U\ <g> V\ acts on Alice's side 
and U 2 (g> V 2 on Bob's. In this sense the product state is 
invariant under the group G <g> H of local unitaries, and 
we can apply the methods developed below to compute 
various entanglement measures for it. 

Computing the commutant (G ® H) 1 is easy, because 
we do not have to look at the Alice— Bob splitting of the 
Hilbert space. In fact, we can invoke the "Commutation 
Theorem" for von Neumann algebras to get 



(G ® H)' = G'®H' 



(22) 



where the notation on the right hand side is the tensor 
product of algebras, i.e., this is the set of all linear combi- 
nations of elements of the form A ® B where A £ G' acts 
on the first two and B £ H' acts on the second two fac- 
tors of Tti®Tt 2 (g)lCi®]C 2 . In particular, if G' and H' are 
abelian, so is G' ® H' and we can readily compute the 
minimal projections, which correspond to the extremal 
invariant states: if p a are the minimal projections of G' 
and qp are those of H', then the minimal projections of 
G' <X> H' are all p a ®qf3- 

Partial transposition also behaves naturally with re- 
spect to tensor products, which implies that (G £g> H)~ = 
G <gi H , and allows us to compute in a simple way the 
(G (8> H )-invariant states with positive partial transpose 
from the corresponding data of G and H . However, for 
the determination of PT> no such shortcut exists. 

We illustrate this in the example, which we will also 
use for the counterexample to additivity of the relative 
entropy of entanglement announced in the Introduction. 
For this we take G — H = UU, with a one-particle space 
Til = Ti 2 = K\ = K,<z = <C d , for any dimension d < 00. 
The extreme points of the state space of G' are given by 
the normalized projections 



1 



P± 



d(d±l) 



(1±F) 



(23) 



Hence the state space of the abelian algebra (G' <g> H') is 
spanned by the four states p Sl (8 p S2 , si,s 2 — ± and is a 
tetrahedron. A convenient coordinate system is given by 
the expectations of the three operators 



F 2 
F\ 2 



F®1 

F®F 



(24) 
(25) 
(26) 



The four extreme points are then on the edges of the unit 
cube: p Sl ® p S2 has expectation triple (si, s 2 , sis 2 ). This 
is drawn in Figure 0. 
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FIG. 4. State space for UUVV-mvaxiant states plotted 
for dimension d = 3. 



pair 2, i.e., 12|34. We claim that upon twirling this pure 
state becomes p#. For this we only need to evaluate the 
expectations of the three operators (|2^), and compare 
with those of p#. Clearly, '5 is a symmetric product 
(Bose-) vector with respect to the total flip F12, hence 
this operator has expectation 1. The expectations of F\ 
and F2 are equal to 

(9\FiV) = T^(ijij\(F® l)\klkt) 
d z £ — ' 

= ^ E m\ekw = ~ ■ 

i,j,k,t 

Since the other four extreme points are separable as ten- 
sor products of separable states, we conclude that all 
ppt-states are separable in this example, so the solid in 
Figure |I] also represents the separable subset. 

Example 8: Tripartite symmetry: U <g> (U ® £7). 



The idea of symmetry can also be used to study multi- 
partite entanglement. A natural choice of symmetry 
group is the group of all unitaries of the form U ® U ® U, 
The resulting five dimensional state space has been stud- 
ied in great detail in [Q. This study also has a bipartite 
chapter, where this group is considered as a group of local 
unitaries U ® (U ® U) in the sense of the present paper. 
The set of separable states is strictly smaller than the set 
of states with positive partial transposes. However, if we 
enlarge the group to include the unitary 1 ® F, the two 
once again coincide, forming a tetrahedron. 



The extreme points are special instances of product 
states: when p, a are t/£/-invariant states with flip expec- 
tations fi and f%, respectively, the product state p<g>a has 
coordinates (/1, $21 fifz)- Hence the manifold of product 
states is embedded in the state space as a piece of hy- 
perboloid. Partial transposition turns the flip operators 
(^J) into their counterparts using JF instead of F. Hence 
the operators with positive partial transposes are rep- 
resented in the diagram by a tetrahedron with vertices 
(0,0,0), (0,d,0), (d,0,0), and {d,d,d 2 ). The intersec- 
tion, i.e., the set of states with trace equal to one and 
positive partial transpose (represented in Figure | as a 
solid) is a polytope with the five extreme points (0, 0, 0), 
(0, 1, 0), (1, 0, 0), (1, 1, 1), and, on the line connecting the 
origin to the point (d, d, d 2 ), the point (1/d, 1/d, 1). The 
density operator corresponding to this last point is 

p# = -^fp+ ®p+ + -^i~p- ® p- ■ ( 2? ) 

It turns out that p# is separable: Let $ = d^ 1 / 2 ^2 k \kk) 
be a maximally entangled vector, and consider a pure 
state with vector ^ = <&Aiice <8> ^Bob- Note that this is a 
tensor product with respect to the splitting Alice— Bob, 
i.e., 13 1 24 rather than the splitting between pair 1 and 



III. ENTANGLEMENT MEASURES AND 
ADDITIVITY 

A. Entanglement of Formation and the convex hull 
construction for functions 

The entanglement of a pure state is well described by 
the von Neumann entropy of its restricted density oper- 
ator. Thus for a pure state p = |4 r )(^ r | such that ^ is 
expressed in Schmidt form as \& = J2k V^k e k <S> e' k , we 
have 

E(P) =&( Cfc ) with ( 28 ) 

k 

r,(t)=-tlog(t). (29) 

The entanglement of formation is a specific extension of 
this function to mixed states. The extension method is 
a general one, known as the convex hull construction for 
functions, and since we will need this construction for 
stating our main result, we will briefly review it. 

So let if be a compact convex set, let M C K be an 
arbitrary subset, and let / : M — > 1R U {+00}. We then 
define a function cof : K — ► H U {+00} by 
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co/(ar) = inf { ]T Ai/(*i)|*< G M, £ A;s, = x} , (30) 

i i 

where the infimum is over all convex combinations with 
Aj > 0, Ai = 1, and by convention the infimum over 
an empty set is +oo. The name "convex hull" of this 
function is due to the property that co/ is the largest 
convex function, which is < / at all points, where / is 
defined. Another way of putting this is to say that the 
"supergraph" of co/, i.e., {(x,r) € K x H|r > co/(x)}, 
is the convex hull (as a subset of K x H) of { (x, r) e 
K x B\x G M, r > f(x)}. 

In this notation, the usual definition || of entangle- 
ment of formation is then 



E F (p) = (coE)(p) 



(31) 



where on the right hand side E is understood as the 
function (|2^) defined only on the submanifold M C 
S(H\ <8> H2) of pure states. 



B. Relative Entropy of Entanglement 

Another measure of entanglement, originally proposed 
in is based on the idea that entanglement should be 
zero for separable density operators (see Equation (||)), 
and should increase as we move away from T>. Such a 
function might be viewed as measuring some kind of dis- 
tance of the state to the set V of separable states. If one 
takes this idea literally, and uses the relative entropy Jn| 



S(p, a) = trp(log p - log a) 



(32) 



to measure the "distance" , one arrives at the relative en- 
tropy of entanglement 



E RE (p) =M{S{p,a)\a eV} . 



(33) 



Initially, other distance functions have also been used to 
define measures of entanglement. However, the one based 
on the relative entropy is the only proposal, which coin- 
cides on pure states with the "canonical" choice described 
in Equation (p8|). Since -Ere is easily shown to be convex, 
it must be smaller than the largest convex function with 
this property, namely Ep. Another reason to prefer rel- 
ative entropy over other distance-like functionals is that 
it has good additivity properties. The hope that -Ere 
might be additive was borne out by initial explorations, 
and has become a folk conjecture in the field. However, 
we will give a counterexample below. 



C. Additivity 

A key problem in the current discussion of entangle- 
ment measures is the question, which of these are "ad- 
ditive" in the following sense: if p, a are bipartite states 



on the Hilbert spaces H\ ® H2 and JC\ (g) IC2, then p(g> a 
is a state on Hi <8> H2 <£> K-i <£> K-2- After sorting the fac- 
tors in this tensor product into spaces Hi , JC\ belonging 
to Alice and 7^2)^2 belonging to Bob, we can consider 
p <g> a as a bipartite state on (Hi <8 ICi) ® (H2 <8> £2). 
This corresponds precisely to the situation of a source 
distributing particles to Alice and Bob, p (£> a, and simi- 
lar larger tensor products, being interpreted as the state 
obtained by letting Alice and Bob collect their respective 
particles. Additivity of an entanglement measure E is 
then the equation 



E(p <g) a) = E(p) + E{a) 



(34) 



We speak of subadditivity if "<" holds instead of equal- 
ity here. Both Ere and Ep are defined as infima, and 
for a product we can insert tensor products of convex 
decompositions or closest separable points into these in- 
fima, and use the additivity properties of entropy to get 
subadditivity in both cases. It is the converse inequality, 
which presents all the difficulties, i.e., the statement that 
in these minimization problems the tensor product solu- 
tions (and not some entangled options) are already the 
best. 

Additivity of an entanglement functional is a strong 
expression of the resource character of entanglement. Ac- 
cording to an additive functional, sharing two particles 
from the same preparing device is exactly "twice as use- 
ful" to Alice and Bob as having just one. Here preparing 
two pairs means preparing independent pairs, expressed 
by the tensor product in (^). It is interesting to in- 
vestigate the influence of correlations and entanglement 
between the different pairs. On the one hand, Alice and 
Bob might not be aware of such correlations, and use 
the pairs as if they were independent. On the other 
hand, they might make use of the exact form of the 
state, including all correlations. Is the second possibility 
always preferable? Entanglement functionals answering 
this question with "yes" have a property stronger than 
additivity, called strong super additivity. It is written as 



E(p) > E{p H ) + E(p K ) , 



(35) 



where p is a density operator for two pairs (four parti- 
cles altogether), and pu and pjc are the restrictions to 
the first and second pair. An entanglement functional 
satisfying this as well as subadditivity is clearly additive. 
Since additivity is already difficult to decide, it is clear 
that strong superadditivity is not known for any of the 
standard measures of entanglement. 

One case of strong superadditivity is satisfied both for 
Ep and Ere, and we establish this property here in or- 
der to get a more focused search for counterexamples 
later on: We claim that ( j§5| ) holds, whenever px, is sep- 
arable, in which case, of course, the second term on the 
right vanishes (as a special case of additivity, when p/c 
is even a product, this was noted recently in |18|). We 
will show this by establishing another property, called 
monotonicity: for both E = Ep and E = Ere, we claim 







E{ P ) > E{ PH ) . (36) 

Monotonicity for £r,e follows readily from a similar 
property of the relative entropy: if pn , a-H denote the 
restrictions of states p, a to the same subsystem, then 
^{PH, <J n) < <t). But if a is separable in (|33"|), then 
so is its restriction au- The infimum over all separable 
states on Hi £g) H 2 is still smaller, hence monotonicity 
holds. 

Monotonicity for Ep is similar: We may do the reduc- 
tion in stages, i.e., first reduce Alice's and then Bob's 
system, and because Ep symmetric with respect to the 
exchange of Alice and Bob, it suffices to consider the case 
of a reduction on only one side, i.e., the restriction from 
Hi ® (H 2 ®K 2 ) to Hi®H 2 . 

Let p be a state on Hi (£> (H 2 ®K. 2 ) and p' its restriction 
to Hi <S>H 2 . 

Consider the states Si on the larger space appearing 
in the minimizing convex decomposition of p, and let s[ 
denote their restrictions to Hi®H 2 . Of course, both Sj 
and s[ have the same restriction to the first factor Hi. 
Hence 

E F (p)=J2^f«), (37) 



will try to reduce the computation of the infimum (|3C|), 
which involves a variation over all convex decompositions 
of x in the high dimensional set K to a computation, 
which can be done entirely in PK. To this end, we define 
the function e : PK ->RU {+00} by 

e(x) = inf {/(s)|s S M, Ps = x} , (38) 

again with the convention that the infimum over the 
empty set is +00. Then the main result of this subsection 
is that, for x £ PK, 

cof(x) — co e(x) , (39) 

where the convex hull on the left is defined by (|30|), but 
the convex hull on the right is now to be computed in 
the convex subset PK. 

We thus arrive at the following recipe for computing 
the entanglement of formation of G-invariant states: 

• Find, for every state p £ PS, the set M p of pure 
states a such that Per = p. 

• Compute 

e(p) := inf {E{a)\a £ M p ) . (40) 



where /(<r) denotes the von Neumann entropy of the re- 
striction of a state a to Hi, and X^^i s i = P ' ■ Because 
the entropy of the restriction is a concave function, the 
value of the sum ( |37| ) can be made smaller by replac- 
ing each s'i with a decomposition into pure states on 
Hi ® H 2 . Minimizing over all such decompositions of 
p' yields Ep(p'), which is hence smaller than Ep(p). 



IV. ENTANGLEMENT OF FORMATION 

A. Simplified computation 

Our method for computing the entanglement of forma- 
tion can also be explained in the gene ral se tting of the 
convex hull construction in Subsection III A , and this is 



perhaps the best way to see the geometrical content. So 
in an addition to a subset M c K of a compact convex set 
and a function / : M — > H U {+00}, consider a compact 
group G of symmetries acting on K by transformations 
ajj : K — ► K, which preserve convex combinations. We 
also assume that ajjM C M, and j{ajjs) — f(s) for 
s £ M. All this is readily verified for au(A) = UAU* 
and / the entanglement defined on the subset M C K of 
pure bipartite states. Our task is to compute cof(x) for 
all G-invariant x £ K, i.e., those with ajj(x) = x for all 
U £ G. 

Since the integral with respect to the Haar measure is 
itself a convex combination, we can define, as before, the 
projection P : K — > K by Px = J dU ajj(x). The set of 
projected points Px will be denoted by PK. Usually, this 
will be a much lower dimensional object than K, so we 



• For later use try to get a good understanding of the 
pure states achieving this minimum. 

• Compute the convex hull of the function ( f4(i| ) . 

The following simplifications are sometimes possible: 
first of all, all pure states in an orbit of G give the 
same value of E, hence we may replace M p by a suit- 
ably parametrized subset containing at least one element 
from every orbit. At this stage it is sometimes already 
possible to discard further states, in favour of others "ob- 
viously" giving a smaller value of E. The final stage is 
sometimes carried out by showing that the function e 
is convex to begin with, but, as we will see, this is not 
always the case. 

The remainder of this subsection is devoted to the 
proof of Equation (^). We will proceed by showing that 
both sides are equal to 

z = inf { Yl A >/( s *)| s > e M > AjPs ' = x } ■ ( 41 ) 

i i 

Indeed, the only difference between ( fill ) and (|30| ) is that 
in ( |4l| ) a weaker condition is demanded on the Si. Hence 
more Si are admissible, and this infimum is smaller, Z < 
cof(x). On the other hand, if Si satisfying the constraint 
for Z are given, inserting the definition of P produces a 
convex combination giving x, namely the combination of 
the states au(si), labeled by the pair (i,U), and weighted 
with Xi JdU . This convex combination is admissi- 
ble for the infimum defining co/, and gives the value 
Ei^iJdU f(au( Si )) = Ei^ifdU f( Si ) = E<Ai/(»i), 
where we have used the invariance property of / and the 
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normalization of the Haar measure. Hence all numbers 
arising in the infimum ( ftl|) also appear in the infimum 
©, which proves that Z < cof(x), hence Z — cof(x). 
In order to prove the equality Z — co s{x) just note 
that in the infimum (^) the constraint is only in terms 
of Psi, whereas the functional to be minimized involves 
f(si). Therefore we can compute the infimum (|4l]) in 
stages, by first fixing all Ps^ and minimizing each f(si) 
under this constraint, which amounts to replacing / by 
e, and then varying over the Psi, which is the infimum 
defining co e. Hence co e{x) = Z = cof(x). 



B. Extending the computation to some 
non-symmetric states 

It is a basic feature of the convex hull that whenever 
the infimum in (^) is found at a non-trivial convex com- 
bination, there is a "flat piece" in the graph of co/, i.e., 
co/ is also known on the convex hull of the minimiz- 
ing Si |Tl[| . The geometrical meaning of this elementary 
observation is immediately clear from low dimensional 
pictures. It is also easy to prove in general: 

Suppose that AjSj = x is a convex decomposi- 
tion of x (with A; > 0) minimizing J2i ^i/( s i); an d let 
x' = X'^i be another convex combination of the same 
points Si . We claim that this convex combination solves 
the minimization problem for co/(x'), i.e., 



co/K) 



(42) 



Indeed, let x' = ■ pjtj be any convex combination with 
tj G M. Then we can find a small number e > such 
that (A, — eA^) > for all i. Hence 



is a convex combination of elements from M represent- 
ing x. But since the decomposition using only the Sj is 
optimal, we have 



From this we immediately get the claimed optimality of 

x ' = YjjXjSj- 

These remarks are especially useful for the case of en- 
tanglement of formation, for any mixed state the opti- 
mizing convex decomposition necessarily involves several 
terms. Hence any computation of an entanglement of 
formation immediately extends to a larger class of states. 
Therefore, it is of great interest not only to get the value 
of the entanglement of formation for a given mixed state, 
but also to find the set of pure states solving the varia- 
tional problem defining E-p 

The symmetric situation studied in this paper is ex- 
treme in this regard: The minimizing sets are always 



complete orbits of the symmetry group. Therefore we 
get a fairly large set of non-symmetric mixed states for 
which the computations below also give the exact value 
otE ¥ . 



C. Results for G = UU 

In this subsection we will apply the general method to 
computing the entanglement of formation for the states 
of Example 1 . 

In the first step we have to determine the set Mf of 
vectors $ £H®H such that = /. In terms of 

the vector components we get 



(43) 



On the other hand, the reduced density operator has 
components pij = J2k ^ik^jk or, in matrix notation, 
p = <f><i>*. Here we may introduce a simplification due 
to U ® U symmetry, by choosing p diagonal. Note, how- 
ever, that we can not choose the restriction to the second 
system, i.e., $ T i> to be diagonal at the same time with- 
out loss of generality In any case, the eigenvalues of p 
become pa = J^k l^ifcj 2 - Hence the pure state entangle- 



(44) 



where n is the entropy function from ( J2E 

For analyzing the variational problem it is useful to 
consider the contributions of each pair of variables 3>y 
and $ji, and of each diagonal element $>u separately. 
The weights of these contributions are 



ment of <&, which by (|28|) is the entropy of p is 



E(i*x*i) = l>(Ei®tt 



hj = \$ij\ 2 + \®ji\ 2 , for Kj 



(45) 
(46) 



The normalized contribution of one such pair or diagonal 
element to / is 



fa = A^. 1 2«(<i.,/i>„) , for ; 3 
fa = 1 , so that 
/ = Aj j fij . 

i<j 



(47) 
(48) 
(49) 



Similarly, we can write the probability distribution 
pn, . . . , pad as a convex combination of probability dis- 
tributions with respective entropies 



Sit , 



(50) 
(51) 



where we have used the abbreviation i?2(p) = ? 7(p)+ 7 7(l — 
p) for the entropy of a two point probability distribution 
(p, 1 — p) . By concavity of the entropy we have 
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(52) 



To find the lower bound on Sy given /y is just another 
instance of the variational problem we are solving, albeit 
with the considerable simplification that only one off- 
diagonal pair of components of $ is non-zero. This leaves 
the following problem: 

Given two complex variables x, y with the 



constraint la;! 2 



W 



minimize s 



mm- 



1, with 2R(xy) = /, 



Since s is monotonely increasing in |ir| 2 form to 1/2, 
this is equivalent to minimizing |x| 2 , given /. The pairs 
{\x\ 2 , /) compatible with the constraints form the convex 
set 



{(A,/)| |/| < 2yJ\(l - A); 0<A<l} . 

From this we get the minimal admissible |.t| 2 = (1 — 
y/l — / 2 )/2 in the above two variable variational prob- 
lem. Hence 



>ea(/«) = ^a(i (l - ^Wi 



(53) 



This function 82 can be shown to be convex by explic- 
itly computing the second derivative and expanding log- 
arithms in a power series. Combining the bounds (52), 
(|5l|), and ( |53] ) with the convexity of £2, we get 



£(I*X*I) > 5>« £ a(/«) 

i<j 

i<j 

= £2 (/-E A ») 



(xl + yF){U ® J7)|12)(12|({7 ® U)*(xl + yF)* 

with fixed x, y, and arbitrary U , have the same E-p. For 
determining these convex combinations we can drop the 
outer factors, and afterwards shift the operators found 
with (xl + yF) G G'. Let 

T = co{(U ® U)\12)(12\(U ® U)*\U unitary} . (55) 

Clearly, every operator in T is a separable density oper- 
ator with flip expectation zero. Conversely, any operator 
p with these properties may be decomposed into pure 
product states \<f> ® ®ip\. These must also have flip 
expectation zero, which means that <f> _L tjj, so that there 
is a unitary U with <f> eg) ip = (U <E) U)\12). Consequently 

Hence in order to determine whether for a given p we 
can compute E-p(p), we transform it to p by the appro- 
priate (xl + yF) -1 , and then test the separability of p. 
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FIG. 5. e-function for UU invariant states. 



Now suppose that / > 0. Then we can choose just 
a single diagonal entry to be non-zero, and find 
i?(|$)($|) = 0, which is clearly the minimum. However, 
if / < the last equation shows that letting any diago- 
nal entry <f>u ^ decreases the argument of £2 further 
in a range where this function is monotonely decreasing. 
Hence the optimum is choosing all = 0, and allowing 
only two non-zero components $^ and $>ji for some i =/= j. 
This concludes the computation of Ep for [///-invariant 
states (see summary below) . 



However, as noted in Section IV B, knowledge of the 



minimizers for e automatically leads to an extension of 
the computation to some non-invariant states. Let x, y 
be a solution of the two variable variational problem with 
/ = tr(pF). Then the minimizing vector is of the form 



x\l2) + y\2l) = (xl + yF)\l2) 



(54) 



All U C3> [/-translates of this vector will do just as well 
and appear in the minimizing decomposition of the UU- 
invariant state. Hence all convex combinations of the 
density operators 



Let us summarize: 

• For the U ®U -invariant state p with 
tr(pF) = f < 0, we have 



1 



(56) 



independently of the dimension d of the underlying 
Hilbert space. When / > 0, the state p is separable, 
and E F (p) = 0. 

• Let p be a (not necessarily invariant) density oper- 
ator with tr(pF) = f and — 1 < / < 0. Then with 
suitably chosen a, (3 € R 



p=(al + (3F)*p{al + (3F) 



(57) 



is a density operator with tv(pF) = 0. Suppose 
that p is separable. Then formula (|56| ) also holds 
for p. 
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D. Results for G = UU 

The computation of the entanglement of formation for 
Example 2 is already known | jl2[ . The minimizing pure 
states are of the form 



(xI + yF)\U) 



(58) 



with real x, y. 

The extension to non-invariant states works in princi- 
ple similar to the UU-c&se, but for d > 2 it is getting a bit 
more complicated, because the e-function is not convex 
anymore. 



values /i, /2- If now the above procedure works for 
these two states, then we have found an optimal 
decomposition for p and can easily compute the 
entanglement of formation. 



E. Results for OO-invariant states 

Here the extension method of Section IV. B turns out 
to do much of the work. The state space, plotted in 
Figure |^, is separated in four regions. The separable 
square and the three triangles A, B, C. 




(F) 

2 _ 3 4 v ' 

FIG. 6. e-function for UU invariant states for d — 2,3, 4. 
The functions are not convex near the right endpoint for 
d > 3. 
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trfjoF) = / > g, we have 

E F (p) = co(h 2 {j) + (1 - 7 ) log(d - 1) 
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(59) 



with 7 = £ (yj +y /(d-l)(d-f)y. For d > 2 

we need also to compute the convex hull. When 
/ < 4, the state p is separable, and E F (p) 



0. 



Let p be a (not necessarily invariant) density oper- 
ator with tr(pF) = /, 1 < / < d and co(e(/)) = 
e(f). Then with suitably chosen a,j3 G H 



{at + /3F)*p{al + f3F) 



(60) 



is a density operator with tr (pF) = 1. Suppose 
that p is separable. Then formula ( |59| ) also holds 
for p. 

• If / satisfies co(e(/)) < e(/) , the convex hull has 
a flat section between f\ < f < / 2 where /i , fa 
are the two end points of the flat piece satisfying 
co(e(/i/ 2 )) = e(/i/2)- We can always find a convex 
decomposition of p in two states with expectation 



FIG. 7. The state space for OO-invariant states seems to 
split naturally in four regions. The separable square PV and 
the three triangles A, B, C. 



In order to apply the extension method to the UU- 
invariant states, we have to see which states can be writ- 
ten as p = (xl+yF)p(xl+yF)* , with pa separable state 
with flip expectation zero. If we take for p any state at 
the left edge of the separable square, it is clear that we 
will get an OO-invariant state again. The explicit com- 
putation shows that with this method we get Ep(p) in 
the full triangle B. Note that by this construction E-p(p) 
depends only on the expectation (F), and not on (F) 
or the dimension d. Employing similarly the extension 
method for UU we find E-p(p) in the triangle A, getting 
a function depending only on (F) and the dimension, but 
not on (F). 



13 



F. Results for Bell-States 

The Bcll-Statcs were one of the first classes for which 
entanglement of formation could be calculated ||. Of 
course, our method reproduces this result, albeit with a 
more economical decomposition. This is a feature shared 
with the Wootters formula flO[| , It is a natural question 
whether the extension method, applied in this basic ex- 
ample, reproduces the Wootters formula. However, it 
turns out that one gets the result only on state man- 
ifolds of lower dimension. We also did not succeed in 
finding another group of local symmetries, which would 
give Wootter's formula in full generality. 

V. RELATIVE ENTROPY OF ENTANGLEMENT 

A. Simplified computation 

Symmetry simplifies the computation of the relative 
entropy of entanglement dramatically: it reduces the 
variation in ( |33"| ) from a variation over all separable states 
a € T> to those, which are also G-invariant. i.e., when 
p — Pp, we have 

E BE (p)=inf{$(p,a)\aePV} . (61) 

The only ingredients of the proof are the convexity of 
a i— ► S(p, ct), the invariance of relative entropy under (lo- 
cal) unitary transformations of both its arguments, and 
that T> is a convex set invariant under local unitaries. In- 
deed, the properties of T> imply that for any a in the full 
variational problem, Per g PI? C T> is also a legitimate 
argument, and the convexity properties of S show that 
this cannot increase S(p, a). Hence the variation may be 
restricted as in (|6l|). We have listed the ingredients of 
the proof so explicitly, because many variations of -Ere 
may be of interest. For example, the "distance" function 
relative entropy can be replaced by a host of other func- 
tions, like norm differences of any kind. The set T>, too, 
may be replaced, for example by the set of ppt-states, 
as suggested by Rains Q , who also made similar use of 
symmetry. 

A second simplification concerns the computation of 
S(p, a) itself, when both arguments are G-invariant. We 
have seen that G-invariant states can be considered as 
states on the commutant algebra G'. Now the relative 
entropy is defined for pairs of states on arbitrary C*- 
algebras [ fl9| , and the form ( |32| ) involving density matri- 
ces is only the special form valid for a full matrix algebra. 
Since P is a conditional expectation onto G', the result 
does not depend jl9| on whether we compute the relative 
entropy via density matrices, or for the corresponding ab- 
stract linear functionals on G'. Without going into the 
details for general algebras G' here, let us see how this 
helps in the case when G' is abelian, as in most of our 
examples. 



Suppose p a , a = 1, . . . , N are the minimal projections 
of G', and denote by uj a — {tvp a )~ l p a the extremal den- 
sity matrices of PS. Then every p G PS has a unique 
representation as a convex combination 

p=y2 P aL ° a = p° > ( 62 ) 

a a 

where the second expression is at the same time the spec- 
tral resolution of p. If we compute the von Neumann 
entropy — tr(plogp) from this, we find a dependence of 
the result not only on the expectations p a — tr (pp a ), 
but also on the multiplicities tr(p a ), as is quite famil- 
iar from statistical mechanics. On the other hand, the 
fact that relative entropy can be defined for states on 
abstract algebras shows that no such dependence can oc- 
cur for relative entropies. Indeed, the terms involving 
logtr(p a ) from p and a cancel, and we get 

S(p, ct) = Pol ( lo g(Pa) - log(o-a)) , (63) 

OL 

where p a and a a are the respective expectations of p a . 

A typical application is the observation that for UU- 
invariant states the expression for the relative entropy 
of entanglement can be written down in terms of the 
tr(pF), independently of the dimension d of the underly- 
ing Hilbert spaces. 

For UU and [/[/-invariant states the sets of separable 
states are just intervals, and the definition of relative en- 
tropy of entanglement requires a minimization over this 
interval. However, due to a general property of the rela- 
tive entropy, the convexity in both arguments, it is clear 
that the minimum is, in fact always obtained at the end- 
point: if p is the state whose entanglement we want to 
calculate, and a is the minimizing separable state, con- 
vexity implies 

S(p, Xa + (1 - X)p) < A S(p, ct) + (1 - A) Sf>, p) 
= XS(p,a) . 

Hence if there were any separable state on the straight 
line segment connecting p and <r, it would give a strictly 
lower minimum, contradicting the minimality of a. 

For UU the boundary separable state has tr(<r_F) = 0, 
i.e., gives equal weight to the minimal projections. We 
have to compute the relative entropy with respect to a 
state with probabilities (1 ± f)/2, i.e., the function 

eRE (/)=log2-s(i±^,^) , (64) 

where we denote by S(pi, . . . ,p„) — —J2kPk^°EPk the 
entropy of a probability vector (p%, . . . ,p n ). This func- 
tion is plotted in Figure g and determines the relative 
entropy of entanglement of [/[/-symmetric states p via 

Eke(p) = e RB (tr(pF)) . (65) 
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FIG. 8. The relative entropy of entanglement for 

UU-invariant states. 



Similarly, the boundary point a of T> for UU invariant 
states is given by ti(aF) = 1. For general / = tc(pF) 
the minimal projections have weights f/d and (1 — f/d). 
Hence for [/[/-symmetric p, we have £"re(p) = for / < 
1, and 



/ 



logd-(l-^)log(d-l)- S[ J -,1 



f 



(66) 



otherwise. For comparison with the results of |T^], note 
that f/d is the so called maximally entangled fraction of 
P- 

Now we look at OO-invariant states. The state space 
and the separable states are drawn in Figure 0. First 
we look at the state with the coordinates (1,3), which 
is a maximal entangled state. The separable states, that 
are minimizing the relative entropy for this state, are the 
states on the whole line connecting the points (0,1) and 
(1, 1). But now we can find the minimizing separable for 
any state in the whole triangle A. We just have to draw 
the straight line connecting the coordinates of a given 
states with the point (1,3). The intersection with the 
border of PT> is then a minimizer for (1,3) and by the 
properties of the relative entropy of entanglement also 
the minimizer for all states on the connecting line. The 
same argumentation works for the edge point (—1,0) and 
the separable border between (0, 0) and (0, 1) giving us 
all minimizers for the triangle B. The whole triangle C 
has the same minimizer, namely (0, 1). 



B. Counterexample to additivity 

To find a counterexample to the additivity of the rela- 
tive entropy of entanglement we use the group introduced 
in Example 7 . We also know that additivity will hold for 
any states where one of the two independently prepared 
states is separable. So in our example we can restrict to 
the area, where both expectation values of F± and F 2 are 
negative. 



For simplicity we increase the group with -FAiicc^^Bob 
leading us to a smaller commutant only spanned by 1 ®> 
1, J* 1 (g> F,l® F + F ® I. As coordinate system we use 
the expectation values of 



F 

F\2 



1 

2 

F ®F 



(I® F + F® 1) 



(67) 
(68) 



The state space is drawn in Figure 




(F) 



FIG. 
= 3. 



State space for [/C/W&-F-invariant states for 



It is just the intersection of the state space of the 
original group (see Figure ^) with the plane given by 
(Fi) = (F 2 ) = (F). The product states, in the sense of 
additivity, are given by the line ((F), (F) 2 ). 

The counterexample we want to look at, is the state 
referring to the coordinates (—1,1), which is given by 
p_ (g) p_ where p_ denotes the normalized projection on 
the antisymmetric subspace of C ® <D d . From Equa- 
tion ( |65| ) we know the relative entropy of entanglement 
for p_ to be log 2 independent from the dimension d. The 
minimizing state was the state with flip expectation value 
equal to zero now denoted as po. So the expected min- 
imizer for the tensor product would be po (g> po located 
on the quadratic product line with the expectation val- 
ues (0, 0). This one gives us the expected value of log 4 
for the relative entropy. Now we calculate the relative 
entropy between p® 2 = p_ ® p_ and p#. 



S(P^,P#) 



tr(p? 2 log p? 2 -p^ 2 log p # ) (69) 
d- 1 



52 i„„ «>2 



= tr(p»" log pT-pT log. 



2d » ? "> <™' 
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Indeed the minimum must be attained on the line con- 
necting po g) po and p# and it can easily be verified, that 
the minimum always is attained on p#. For d = 2 the 
whole line gives the same value and although there exits 
minimizer not belonging to the product space, additiv- 
ity holds. For d > 2 the expectation values of state p# 
given by (^,1) shift near to the F12 axis and from a ge- 
ometrical point of view closer to p_ <g> p_ . Although the 
relative entropy is not a real kind of geometrical mea- 
sure this intuition did not fail. In these cases the ad- 
ditivity is violated with an amount of log 2 ^J" 1 ^ . For 
very high dimension d we get the really surprising result 
E RE (p- <8> p-) ~ Ejie (P— ) ■ 



VI. CONCLUDING REMARKS 

We have concentrated on just two basic entanglement 
measures. Clearly, there are many more, and for many 
of them the computation can be simplified for symmet- 
ric states. Among these measures of entanglement are 
the "best separable approximation" of a state (2(J, the 
trace norm of the partial transpose |2^]j_the base norm 
associated with V (called cross norm in [E2| and absolute 
robustness in (23)). For distillible entanglement we refer 
to the recent paper of Rains || . Similarly, there is a lot 
of work left to be done carrying out the programme out- 
lined in this pap er fo r all the groups of local symmetries 
listed in Section II A . 
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